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ABSTRACT 

We use a series of high-resolution N-body simulations of the concordance cosmology to 
investigate the redshift evolution since z = 1 of the properties and alignment with the 
large-scale structure of haloes in clusters, filaments, sheets and voids. We find that: 
(i) Once a rescaling of the halo mass with M*(z), the typical mass scale collapsing 
at redshift z, is performed, there is no further significant redshift dependence in the 
halo properties; (m) The environment influences the halo shape and formation time 
at all investigated redshifts for haloes with masses M < M*; and (iii) There is a 
significant alignment of both spin and shape of haloes with filaments and sheets. In 
detail, at all redshifts up to z = 1: a) Haloes with masses below ~ M, tend to be more 
oblate when located in clusters than in the other environments; this trend is reversed 
at higher masses: above about M^, halos in clusters are typically more prolate than 
similar massive halos in sheets, filaments and voids, b) The haloes with M > in 
filaments spin more rapidly than similar mass haloes in clusters; haloes in voids have 
the lowest median spin parameters; c) Haloes with M < tend to be younger in 
voids and older in clusters; d) In sheets, halo spin vectors tend to lie preferentially 
within the sheet plane independent of halo mass; in filaments, instead, haloes with 
AI < tend to spin parallel to the filament and higher mass haloes perpendicular 
to it. For halo masses M > A/*, the major axis of haloes in filaments and sheets is 
strongly aligned with the host filament or the sheet plane, respectively. Such halo-LSS 
alignments may be of importance in weak lensing analyses of cosmic shear. A question 
that is opened by our study is why, in the < 2 < 1 redshift regime that we have 
investigated, the mass scale for gravitational collapse, M*, sets roughly the threshold 
below which the large-scale structure environment either begins to affect, or reverses, 
fundamental properties of dark matter haloes. 

Key words: cosmology: theory, dark matter, large-scale structure of Universe - 
galaxies: haloes - methods: N-body simulations 



1 INTRODUCTION 

Numerical simulations of concordance cosmology have 
shown that properties of dark matter haloes do not depend 
only on the mass of the halo, as suggested by prior ana- 
lytical work based on the excursion set theory for struc - 
ture evolution (e.g. iBond et al.lll99ll : [Lacev fc CoQll993h . 
Rather, they also d epend on the environment in which 



the halo resides (e . g. Gao et al 



iHarker et al.ll2006l : IWang et al 



■20051: IWechsler et al.ll2006l : 
20071 : iHahn et al.ll2007l . Pa- 



per I in the following). This dependence on environment 
is quite significant at redshift zero for low mass haloes, 
typically with masses < 5 x W^^ h~^MQ. In detail, us- 
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ing marked statistics, ISheth fc Tormeiil (|2004l ) find evi- 
dence that haloes of a given mass form earlier in dense re- 
gions. High resolution simulations confirm this finding and 
quantify it as a function of halo mas s llGao et al 



Croton et all 120071: IHarker et all I2OO6I : iReed et al 



2005 



2007 



Maulbetsch et all bOOTf T At the same time it has become 



clear that also other halo properties as concentrati on and 
spin correlate with l ocal e nvironment llAv ila-Rcc se"et al. | 



2005; Wcchslcr ct al. 2006|; iBett et all 120071: iMaccio et al 



2006: Wetzel et al., .20071 : iHahn et al]|2007^ . iGao fc Whiti 



(2003) ind that haloes with, e.g., high spin-parameter or 
formation time, tend to be more strongly clustered than 
younger and low-spin haloes. It is possible that this environ- 
mental dependence of halo properties has also an impact on 
the baryonic galaxies. Galaxy properties in the local Uni- 
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verse are known to vary systematically with environment 
(e.R. [Presslerl Il980l : iKauffmann et aL I I2OO4I : iBlanton et alj 
I2OO5I ). 

In this paper, we extend the investigation of the prop- 
erties of dark matter haloes as a function of environment 
to high redshifts. In particular, we follow the definition of 
environment that we presented in Paper I, which associates 
haloes to four classes with different dynamical properties: 
voids, sheets, filaments and clusters. These four environ- 
ments are identified on the basis of a tidal stability criterion 
for test p articles which is inspired by the Zel'dovich approx- 
imation (Z erdovichlll97d ). We find that, at each redshift, all 
investigated properties of haloes show some correlation with 
mass and environment, and that the redshift dependence of 
halo properties with mass is removed when such properties 
are investigated as a function of the rescaled mass M/M*, 
where M* is the typical mass-scale collapsing at each epoch. 

The paper is organised as follows. In Section [J] we 
briefiy summarise the specifics of our N-body simulations, 
the definitions of the halo properties that we study and the 
definitions of the four environments - clusters, sheets, fila- 
ments and voids. We present the results in Section |3] and 
summarise our conclusions in Section |¥1 



2 NUMERICAL SIMULATIONS AND 
DEFINITIONS 

We use the three high-resolution cosmological N-body sim- 
ulations described in more detail in Paper I, which wer e 
obtained with the tree-PM code Gadget-2 (|Springelll2005h . 
These simulations are used to follow the formation and evo- 
lution of large-scale structure in a flat ACDM cosmology 
with matter density parameter f2m = 0.25, baryonic contri- 
bution S7b = 0.045 and a present-day value of the Hubble 
constant Hq = 100 ftkms"^ Mpc~^ with h — 0.73 with an 
initial power spectrum normalised to have erg = 0.9. Each 
simulation consists of 512^ coUisionless dark matter particles 
in periodic boxes of sizes Li = 45 Mpc, L2 = 90 Mpc 
and L3 = 180 Mpc, respectively. The corresponding par- 
ticle masses are 4.7 x lO'', 3.8 x 10® and 3.0 x 10^ /i'^Mq for 
the three boxe s. Initial condition s were generated using the 
Grafic2 tool (|Bertschingeill200ll 'l . Particle positions and ve- 
locities were saved for 30 time-steps, logarithmically spaced 
in expansion parameter a between z = IQ and z = Q. The 
mass-range of these three simulations allows us to resolve 
haloes with masses below M» up to redshifts 2 < 1. 

The halo properties that we investigate are formation 
redshift, shape, and spin parameter. The formation redshift 
and shape parameters are defined as in Paper I; we adopt 
however a slightly different approach to measure the halo 
spin parameter than in our previous work. We summarise 
our definitions below. 

2.1 Halo Catalogues 

Haloes were identified in each sn apshot using the standard 
friends-of-friends (FOF, |d avis et al.. .1985, ) algorithm with 
a linking length equal to 0.2 times the mean inter-particle 
distance. Haloes that are well-resolved in each of the three 
simulations are then combined into one single catalogue. Un- 
relaxed systems were identified and deleted from the halo 



catalogues. These unrelaxed systems are mainly contributed 
by close-pair halo configurations which are spurio usly linked 
into one single halo. To identify them, we follow iBett et al] 
(2007) and define the virialization parameter 



where K is the total kinetic energy including the Hubble 
fiow with respect to the centre of mass and V is the total 
potential energy of the isolated FOF halo. The potential is 
computed using a tree for groups with more than 5000 par- 
ticles, and via direct summation for smaller halos. The virial 
theorem states that the time average of a vanishes for any 
isolated relaxed object. However, in falling material exerts 
a sur face pressure such that a < (|Hetznecker fc BurkertI 
I2OO6I ). In addition, structures that are gravitationally bound 
have a > —3. In order to exclude accidentally linked unviri- 
alised structures or haloes that are just about to merge, it 
suffices to fix bou nds on a. In order to directly compare with 
iBett et al.ll|2007l ). we make the same choice of |q:| < 1/2 that 
was adopted by those authors to set the threshold between 
virialized and non-virialized structures. 

Finally, we also exclude from our halo catalogues all 
those structures where the distance between the centre of 
mass rcM and the most bound particle tmb exceeds a fixed 
fraction / = 0.25 of the largest distance between a particle 
in the halo and the centre of mass Tmax, i.e. / = |rcM — 

rMB|/|rniax|. 

The cleaning of the halo catalogues has a strong effect 
on the spin parameter distribution but only a minor influ- 
ence on the other quantities that we study in this paper. 

2.2 Formation Redshift 

For each halo at redshift z, we identify a progenitor aX Zp > z 
by identifying particles that are contained in both haloes. 
The main progenitor is then chosen to be the most massive 
halo at each redshift that contributes at least 50 per cent of 
its particles to the flnal halo. We then deflne the formation 
redshift Zform to be the epoch at which a main progenitor 
which has at least half of the flnal mass flrst appears in the 
simulation; speciflcally, Zform is found by linearly interpolat- 
ing between simulation snapshots in log z to find the point 
where half of the given halo mass is accumulated. 

2.3 Halo Shape 

In order to determine the shape of haloes, we use the mo- 
ment of inertia tensor 

Ijk = m'^[rfSjk — XijXi^k) , (2) 

i 

where m is the particle mass, = \ {xi,i, Xi,2, Xi^3)\ is the dis- 
tance of the i-th particle from the centre of mass of the halo 
and Sjk denotes the Kronecker symbol. Given the lengths of 
the principal axes of inertia h ^ h h, we then use the 
following definitions of sphericity S and triaxiality T: 

I- P — P 

We find that a minimum of 500 particles per halo guarantees 
numerically reliable estimates of the shape parameters. 
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2.4 Halo Spin Parameter 



We estimate tlie spin pa rameter llPeebled 19691 ') of a iialo 
usme tlie simplified form (jBullock et al.ll200ll l 



A' 



I J vir 



(4) 



Here all quantities with the subscript "vir" (angular mo- 
mentum, mass and circular velocity) are computed within a 
sphere of radius R^ir around the most bound particle enclos- 
ing a mean density of A{z)pc{z), where Pc{z) is the critical 
density, and A{z) the density parameter according to the 
spherical collap se model. This density p arameter can be ap- 
proximated by l|Brvan fc Normanlll998l ): 



A{z) = 18% +82 f{z)~ 39 f{z), 



with 



(1 + 2)3 -f- 



(5) 



(6) 



in a fiat cosmology. Under the assumption that the halo 
is in dynamical equilibrium, V^i,. = GMvir/ Rvir, the spin 
parameter can be rewritten as 



A' 



3/2 ■ 



(7) 



Systematic numerical artefacts were found to be negligible 
for haloes consisting of at least 300 particles. The clean- 
ing of the halo catalogues, as described in Section [2.11 has 
a strong influence on the distribution of A'. The probability 
for two haloes of similar mass to be erroneously linked by the 
halo finder grows both with the environmental density and 
decreasing mass of the haloes. This leads to an increasing 
component of unrelaxed structures of low mass in the spin 
distribution for which the virialisation conditions are not ful- 
filled. Furthermore, the angular momentum J is dominated 
by the orbital angular momentum of the pair rather than 
the intrinsic spin of either one of them. Exclusion of unre- 
laxed objects removes the tail of these apparent high-spin 
haloes with A' > 0.1. 



2.5 Environmental Classification 

We employ the definition of environment that was described 
in detail in Paper I. This classification is based on a local 
stability criterion for the orbits of test particles in the mat- 
ter distribution at a fixed epoch. A series expansion of the 
equation of motion for a test particle in the smoothed matter 
distribution gives a zero order acceleration and a first order 
deformation term that is determined by the tidal field ten- 
sor, i.e. the Hessian of the peculiar gravitational potential. 
We then classify the contraction/expansion behaviour of this 
first order term by the number of its positi ve/negative eigen - 
values. In analogy with Zel'dovich theory llZel'dovichll 1970f ) . 
the application of this criterion provides four environmen- 
tal categories defined by the signs of the three eigenvalues 
Ai < A2 ^ A3, namely: 



(i) clusters (Ai,2,3 ^ 0), 

(ii) filaments {X2.S 0, Ai < 0), 

(iii) sheets (A3 0, Ai,2 < 0), and 

(iv) voids (Ai,2,3 < 0). 



This definition relies on one free parameter, the length scale 
Rs of the Gaussian filter that is used to smooth the matter 
distribution before obtaining the gravitational potential on 
a grid. As done in Paper I, we fix the smoothing scale at 
redshift z = to he Rs — 2.1 h^^ Mpc, as this value pro- 
vides the best agreement between the outcome of the orbit- 
stability criterion and our a-posteriori visual classification of 
the different environments. The smoothing length scale Rs 
is related to the mass Ms contained in the Gaussian filter 
at mean density p via Ms = {2%)'^^^ pRs', thus, a smoothing 
scale Rs = 2.1 Mpc corresponds to Ms = lO" Mq. 

A natural mass scale at any given redshift is given by the 
typical mass scale for collapse Mt{z) defined as follows. A 
spherical top-hat perturbation collapses when its linear over- 
density exceeds a value of 5c — 1.686 with only a weak de - 
pendence on cosmological parameters (e.g. lEke et al.|[l996l ). 
The variance of linear density fiuctuations at a given mass 
scale M is related to the linear power spectrum P{k, z) at 
redshift z by 



dk k'^ P(k,z)WTuik,M), 



(8) 



where WTH{k, M) is the Fourier transform of a spherical top- 
hat window function of comoving size R = (3M / 4npy^^, 
and p is the comoving mean mass density of the universe. 
At a given redshift, the typical mass scale M, (2) to collapse 
from a la fluctuation is hence given by the implicit solution 
of 

a{M,,z)=Sc. (9) 

The mass scale M*(2) at 2 = is 5.5 x 10^^ /I'^M©. 

When computing the environmental classification at 
redshifts 2 > 0, there are two possible approaches that can 
be followed: i) To keep the smoothing scale Rs (Ms) fixed 
to the 2 = value: the environment is thus defined over 
typical separations of a few Mpc in comoving space; or li) 
To vary the smoothing scale. In particular, a natural choice 
for a redshift-dependent smoothing scale is to maintain the 
ratio Ms/M, fixed for the Gaussian filter. The respective 
mass scales Mt{z) for the high-2 simulation snapshots in- 
vestigated in this paper, i.e., 2 — 0.49 and 2 = 1.05, are 
1.2 X 10^^ and 1.9 x 10" h'^ Mq. Fixing the ratio Ms/M* 
maintains the smoothing on scales of order ~ 2M, at all 
redshifts. 

The resulting classifications for both AIs — const, and 
Ms/Mf = const, at redshifts 2 = 0, 0.49 and 1.05 are shown 
in Figure [T] using different colours for the cluster, sheet, 
filament and void environments. We observe some differ- 
ences between the two smoothing approaches. With a fixed 
smoothing length Rs = const., shown in panels b) and d), 
the regions classified as voids, sheets and filaments remain 
virtually unchanged since 2 = 1, while the cluster environ- 
ments grow substantially in size, from 2 = to higher red- 
shifts, and extend well into the filaments at 2 = 1. With 
the constant Ms/M* smoothing, shown in panels c) and 
e), a much larger number of individual haloes change envi- 
ronment with redshift: at the resolution of our simulations, 
very few haloes are detected at 2; = 1 in void regions, while 
many haloes are associated at the same redshift to rela- 
tively small cluster environments. In Tables [T] (for fixed Ms) 
and [2] (for fixed Ms/M,) we show the fraction of haloes at 
2 = that change their environmental class from 2 = 1.05 
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Figure 1. Halo environment classification as a function of redshift in projected slices of 15 /i^^ Mpc for tfie 90 h~^Mpc box. Only 
for presentation purposes, i.e., to increase the contrast in the figure, we represent with a dot haloes down to structures with only 10 
particles. The four environments are identified by colour: voids (orange), sheets (green), filaments (blue) and clusters (red). Panel a) is 
the slice at 2 = 0; the smoothing mass scale is Ms = 10^'' Mq. Panels b) and c) are snapshots at 2 = 0.49; panels c) and d) are at 
2 = 1.05. In panels b) and d) the smoothing scale Ms adopted for the environmental classification is as at 2 = 0; in panels c) and e) it 
is the ratio Ms/M« that is kept fixed. 



through z — 0.49 to z = by following the main progenitors 
of each halo with a minimum mass of 10^^ Mq in the 
box. For a fixed smoothing mass scale, indeed 
less haloes change their environmental class as the density 
contrasts between the environments grow through gravita- 
tional collapse. 

The fraction of volume attributed to each of the four 
environments as a function of redshift is shown for both 
smoothing approaches in Figure [H For a fixed smoothing 
mass Ms = const., the density field asymptotically ap- 
proaches Gaussianity with increasing redshift and thus the 
expected volume fr actions for the four environments (cf. 
lDoroshkevichlll970l ). The behaviour is very different with 



the constant Ms/M« smoothing. The volume occupied by 
the unstable structures (sheets and voids) increases with 
redshift, while the volume fractions of the stable structures 
(clusters and filaments) decrease compared to the values at 
z = 0. 

The redshift evolution of the median value of the 
smoothed overdensity field as measured at the centres of 
haloes in the four environments is shown in Figure O For 
a constant smoothing scale, the median overdensities grow 
faster than expected in linear perturbation theory as \5\ > 1 
in clusters and voids. With the Ms/M^ = const, smoothing, 
however, the median overdensities are found to be essentially 
constant at all redshifts. In both smoothing approaches, the 
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Table 1. Environmental transition matrix for the main progen- 
itor branch of haloes with masses M{z = 0) > 10" /i-i Mq be- 
tween z = and z = 0.49 (upper half) and z = 1.05 (lower half). 
Matrix entries are given in per cent of the total number of haloes 
at z = 0. Non-diagonal elements represent haloes that change 
classification. Environment is determined with Mg = const, (i.e., 
Rs = 2.1 Mpc at all redshifts). 
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Figure 2. Volume fractions for the 90 fc~^Mpc box over red- 
shift. Thick black lines represent the fractions obtained using a 
smoothing scale constant with redshift, and grey lines indicate the 
corresponding fractions obtained when keeping the ratio M/A/» 
constant. Thin dotted black lines represent the values predicted 
for a Gaussian field (42 per cent for sheets and filaments, and 8 
per cent for voids and clusters). 
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Figure 3. Median overdensity in cluster (dot-dashed), filament 
(long dashed) , sheet (short dashed) and void (solid) environments 
for a constant smoothing scale Ms (top) and for a constant ratio 
Ms/M, (bottom) with redshift. The smoothing scale at redshift 
zero is Mg = lO'^^ /j-i Mg. 



median overdensity in filaments is constant. The median 
overdensities, smootlied on scales of Rs = 2.1 Mpc at 
redshift zero, are 5 = —0.79, —0.55, 0.28 and 4.44 in voids, 
sheets, filaments and clusters, respectively. 

We note that with both smoothing approaches, and 
most relevantly when adopting a constant Ms/M* ratio for 
the smoothing, the haloes in the immediate surroundings 
of the clusters a,t z — 1 are classified as filaments/sheets 
at this redshift, but they make the transition to the cluster 
environment by 2 = 0. This allows us to rigorously iden- 
tify, and thus study the properties of, the haloes in these 
intermediate-density "transition regions", before their ulti- 
mate migration into the deeper potential wells of rich clus- 
ters at z = 0. We plan to investigate these haloes in future 
work. 



3 THE REDSHIFT EVOLUTION OF HALO 
PROPERTIES IN DIFFERENT 
ENVIRONMENTS 

3.1 Mass Functions 

The choice of smoothing scale with redshift has an impact on 
the analysis of the redshift evolution of the halo properties 
in the different environments. Starting with the halo mass 
functions, shown in Figure |4] for the cluster, sheet, filament 
and void environments at 2 = 0, 0.49 and 1.05, there is a sub- 
stantial change in their shapes when using one or the other 
of the smoothing approaches. Adopting a constant Ms/M, 
for the smoothing scale leads to a substantial increase in 
low mass haloes that are classified to be in clusters relative 
to the other environments. In voids, sheets and high mass 
filaments, the mass functions are higher when smoothing 
with Ms — const, than when adopting a constant Ms/M, 
ratio; the trend reverses for haloes in low mass filaments and 
clusters, for which the mass functions are instead enhanced 
when using the Ms/M« = const, smoothing scale. The in- 
fiexion point on scales of ~ 10^^ ft"^MQ (~ 10^^'^ /i-IMq) 
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Figure 4. Mass function of tlie haloes residing in voids, slieets, 
filaments and clusters at redshifts 2 = (red), z = 0.49 (green) 
and z = f.05 (blue). Abundances in the whole box have been 
rescaled by the corresponding volume fractions. The dotted curves 
are obtained -with Ms/M* = const., the dashed curves with Mg = 
const at 2 > 0.49 and 1.05. The smoothing scale at redshift zero 
is Ms = 10^^ Mq; here the two smoothings coincide and are 
represented by the solid line. 



in the mass function of z = 1.05 {z — 0.49) clusters for the 
Ms/Mt = const, smoothing reflects the increasing abun- 
dance with redshift of isolated small clusters that we men- 
tion in Section [2.51 At all redshifts of our study, the clusters 
have the highest mean number density of haloes, followed 
by fliaments, sheets and voids, respectively. 



3.2 Halo Formation Redshift 

A closer inspection of the dependence of halo properties 
on the adopted smoothing scale shows that, with both ap- 
proaches, there is always a threshold mass-scale Mt below 
which the environmental influence on halo properties either 
switches sign or transitions from being negligible to becom- 
ing increasingly more substantial with decreasing halo mass, 
down to the ~ 10^" /i~^Mq scales which can be probed 
with our simulations. Figure [5] shows, for both smoothing 
approaches, the dependence on environment and mass (in 
units of M, ) of the median halo formation redshift z^^d- 
Errorbars are estimates of the error in the median, and are 
computed as: 



Ax ■ 



So. 84 — Xo.16 



(10) 



where Nh is the number of haloes used to sample the dis- 
tribution of X, and a::o.84 and a;o.i6 denote the 84th and 
16th percentile of the distribution. These values would corre- 
spond to la errors if the underlying distribution were Gaus- 
sian. The la confidence region of the median formation red- 
shift determined from all environments is represented by the 
shaded regions in Figure [5] As discussed in Paper I, this 
overall behaviour, not split by environment, is well fit by a 
logarithmic relation over five decades in mass at z = 0. This 
relation, refiecting hierarchical assembly is also present at 
higher redshifts. We fit a model of the form 



A least-squares fit to all haloes extracted from our three 
simulation boxes provides the following parameters at the 
three redshifts: 



Cl 


= 1.07 ±0.01, 


C2 


= 0.30 ±0.01; 


Cl 


= 1.70 ±0.01, 


C2 


= 0.33 ±0.01; 


Cl 


= 2.47 ±0.01, 


C2 


= 0.34 ±0.01; 



z = 0, 



0.49, 



z = 1.05. 



We note that, in all plots and thus at all redshifts and for 
both smoothing approaches, there is indeed a mass scale Mt 
at which the curves for the four environments meet, indi- 
cating the vanishing of significant environmental influence 
above this mass threshold. Specifically, below Mt, haloes 
form earlier in clusters than in filaments, followed by sheets 
and then voids. This difference in formation time increases 
with decreasing mass below the threshold Mt. 

It is clear from Figure [S] that, in the case of a smooth- 
ing scale that remains constant with redshift, the threshold 
Alt coincides with the mass-scale for gravitational collapse 
M* at z = 0, but strongly departs (and increases relative 
to M*) at higher redshifts. Interestingly, however, in the 
case of the Ms/M, = const, smoothing scale, the thresh- 
old mass is easily identified to lie very close to M, at all 
redshifts. This difference is simply due to the different asso- 
ciation of haloes to the cluster and filament environments in 
the two smoothing approaches that is also observed in the 
mass functions (cf. Figure |4)). Furthermore, we note that 
the z — 0.49 and z = 1.05 relations are very similar to 
the one at z = 0, for which we had already provided the 
analytic fits with environment-dependent slopes for masses 
M < M, in Paper I. The environmental dependence of the 
halo formation redshifts below the M, mass scale that w e 
have found agrees with the analysis of lGao fc White! l|2007h . 
who find that haloes with higher formation redshifts are 
more strongly clustered. 

Interestingly, at the highest redshift of our study {z = 
1.05), and for the Ms/AI, = const smoothing, our simula- 
tions show that haloes with masses M > M* in filaments 
tend to have earlier formation times than haloes of similar 
masses in the cluster environment, i.e., an opposite trend 
than the one observed at all redshifts below M = M* . 

The fact that, in the case of a fixed Ms/M^ smooth- 
ing ratio, the environmental dependence is explicitly seen to 
appear around A4/M, = 1, hints at a physical relevance of 
this scale in establishing the onset of the environmental de- 
pendence of halo properties at all redshifts. This motivates 
us to identify the constant Ms/M* ratio as the more fun- 
damental smoothing scale in our analysis, and thus to use 
this smoothing scale in the remainder of our study of the 
redshift evolution of halo spins, shapes and alignments as a 
function of environment. 



3.3 Halo Spin 

Figure [S] shows the median spin parameter X[^^c^ as a func- 
tion of normalised mass M/M* for haloes at z = 0, 0.49 and 
1.05 in the cluster, sheet, filament and void environments. 
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Figure 5. Median formation redshift in voids, sheets, filaments and clusters as a function of halo mass in units of M,(z) at redshifts 
2 = 0, 0.49 and 1.05. The left panel shows the results when smoothing with a constant Ms; the right panel shows the case of the 
smoothing with a constant Ms/Mt ratio. Errorbars indicate the uncertainty on the median. The shaded area indicates the Icr confidence 
area for the median of the whole sample not split by environment. 




M / 

Figure 6. Median spin parameter in filaments and clusters as a 
function of halo mass in units of M* at redshifts z = 0, 2 = 0.49 
and z = 1.05. Errorbars indicate the uncertainty on the median. 
The shaded grey area indicates the la confidence area of the 
median of the whole sample not split by environment. 



We also plot the results integrated over all environments 
(shaded region in Figure [6)l , for co mparison with previou s 
studies. In agreement with these (e.g. lVitvitska et al.ll2003 ). 
we do not observe any significant evolution of the global spin 
parameter with redshift. Our results also confirm a weak 
mass depe ndence of t h e med ian spin parameter as found at 
z = by iBett et al.l (|2007l ). and furthermore extend this 
result to significantly lower masses. Integrated over all en- 
vironments, the spin-mass relationship is weakly increasing 
up to « lOM, (z), while the highest masses at each redshift 
have again a slightly smaller median spin parameter. 

We find some dependence of the halo spin on environ- 
ment at all redshifts of our analysis. More specifically, high- 
mass (M > M«) haloes in filaments have a higher median 
spin at all redshifts than comparable-mass haloes in the clus- 



ter environments; in the voids, haloes with masses substan- 
tially below Mt{z) spin systematically more slowly (i.e. me- 
dian A' < 0.03) than haloes of similar masses in any other 
environment. We have tested that this trend of void haloes 
having lower median spin parameters persists and actually 
increases when a larger scale is adopted for the smoothing, 
to optimize the identification of the void regions (cf. Paper 
I). At the lowest masses we do not observe any significant 
diff erence between haloes in clusters, filaments or sheets. 

iGao fc Whitl (|2007| ) report that, in their simulations, 
the most rapidly spinning haloes are more clustered than the 
slowest spinning haloes, which is in agreement with our ear- 
lier findings ai z — (Paper I). Our results, after removing 
unrelaxed haloes as described in Section [211] do not support 
a very strong correlation between environment and spin at 
low masses. Still, it is slightly more likely to find the most 
rapidly spinning objects in environments of higher median 
density. 

3.4 Halo Shape 

In Figure [T] we show the median sphericity and triaxiality 
of haloes in the four environments as a function of their 
mass for redshifts z — 0, 0.49 and 1.05. Independent of red- 
shift and environment, haloes tend to be more spherical with 
decreasing mass. Over a large range of masses, haloes at 
z = 1.05 are h owever less spheri c al tha n ha loes of similar 
mass at z = 0; lAvila-Reese et all (|2005l ) and lAUgood et al.l 
(|2006i) find a similar result in their analyses. The mass-shape 
relations, expressed as scale-free functions of M/M*, show 
no significant evolution with redshift up to z ~ 1; the fact 
that the entire redshift evolution of the shapes of haloes is 
driven by the evolution of the mass scale for gravitational 
collapse, A'/«, is also supported by the independent studies 

quoted above. 

Similarly to the z = case iBett et al.ll2007l . Paper I), 
the mass-sphericity and the mass-triaxiality relation of the 
global (i.e., not split for environment) sample follow a bro- 
ken logarithmic relation also at high redshifts. The change 
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Figure 7. Median iialo spliericity (left) and triaxiality (rigtit) as a function of halo mass in units of M* for haloes in the four environments 
at redshift 2 = 0, 0.49 and 1.05. Environment is defined adopting a constant Ms/M, ratio. Errorbars represent the error in the median. 
The shaded grey area indicates the Icr confidence area of the median of the whole sample not split by environment. 



in slope in these relationships occurs around M « M* . We 
detect however a relatively small but systematic difference in 
the sphericity and triaxiality of M < M, halos in low- and 
high-density environments, i.e., the slope in the M < A/* 
regime weakly depends on the environment. Specifically, the 
median sphericity of M < Mi,(z) haloes decreases systemat- 
ically from the cluster environments, to the filaments, sheets 
and voids. At all redshifts, a weak trend is observed for 
haloes with masses below M* to be more oblate in clusters 
than in filaments, and more prolate in sheets than in fila- 
ments; for haloes above M», there is a stronger evidence for 
haloes in filaments to be systematically more oblate than in 
clusters. These environmental differences at low masses are 
observed to be already in place aX z = 0.49; the resolution 
of our simulations is not adequate to properly investigate 
these effects at z — 1.05 (minimum halo mass ~ O.IM,). 



3.5 Halo Alignments 

3.5.1 Halo-LSS alignment 

Extending the analysis of Paper I to redshifts above zero, 
we use the directional information derived from the eigen- 
structure of the tidal field tensor to estimate the alignment 
of halo spin with the dynamical properties of the surround- 
ing environment. Filaments and sheets have a preferred di- 
rection given by the eigenvector corresponding the single 
positive or negative eigenvalue. The eigenvectors indicating 
the direction of the filament as determined from the tidal 
tensor are shown in Figure [S] Given these unit eigenvectors 
V, we compute the alignment angle cos^ = J • v. Figure [9] 
shows the median alignment as a function of mass at red- 
shifts 2 = 0, 0.49 and 1.05. At all redshifts, there is a strong 
tendency for sheet haloes to have a spin vector preferen- 
tially parallel to the sheet, i.e. orthogonal to the normal 
vector. At redshifts up to 0.49, where the errorbars of our 
measurements allow us to investigate trends with halo mass, 
this alignment increases with increasing mass. For filament 
haloes, there is a clear trend with halo mass: (i) haloes with 
masses smaller than about 0.1 Af, have spins more likely 




h 'Mp( 



Figure 8. Unit eigenvectors indicating the direction of the fil- 
aments are shown in black for filament haloes in a slice of 8 
/i~^Mpc in the 180 Mpc box at z=0. The grey symbols 
indicate halo positions regardless of their environment. The di- 
rectional information of these vectors is used to determine the 
alignment of halo spins with the large-scale structure. 

aligned with the filament in which they reside; (ii) haloes 
in the range M « 0.1 Af* to lAf« appear to be randomly 
aligned with respect to the large-scale structure; and (Hi) 
For M > M*, the trend appears to reverse, and more mas- 
sive haloes have a weak tendency to spin orthogonally to the 
direction of the filament at lower redshiftfl 

^ The tendency for haloes above M» to spin orthogonal to the 
host filament, shown in Figure[9]for the Ms/M« = const, smooth- 
ing case, is enhanced when the Ms = const, smoothing is adopted. 
The smoothing scale not only determines the environmental split 
of the halo population, it also affects the scale on which the eigen- 
vectors of the tidal field are computed. When the smoothing is 
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Figure 9. Median alignment angles between the halo angular 
momentum vectors and the eigenvectors pointing in the direction 
of filaments and normal to the sheets, respectively. Difi'erent red- 
shifts are indicated with the three colours. Errorbars indicate the 
error in the median. The dotted line indicates the expectation 
value for a random signal. 



To further explore possible connections between the 
alignment of the large-scale structure and the intrinsic align- 
ment of haloes in the different environments, we search for 
a correlation signal between the LSS and the axis vectors of 
the moment of inertia ellipsoid of the haloes. In particular, 
we use the major axis vector li to define the alignment an- 
gle cos 8 = li • V, where v is again the eigenvector normal 
to a sheet or parallel to a filament. The resulting median 
correlation is shown in Figure [TO] We find no alignment for 
halo masses M < O.lAf*; however, in both the filaments 
and the sheets, the halo major axis appears to be strongly 
aligned with the LSS for masses above about a tenth of A'h . 
The strength of the alignment grows with increasing mass. 
This is possibly to be expected, especially for the most mas- 
sive haloes, since their shape might infiuence the potential 
from which the eigenvectors are derived. Adopting a fixed 
smoothing scale Ms results merely in a shift of the relations 
shown in Figure 1101 

Results similar to ours concerning the alignments of 
shapes and spins with the LSS, and the transition of align- 
ment orientation at M, i n the filaments, are reported by 
lAragon-Calvo et al.l (|2007l ) for z — haloes using a defini- 
tion of environment that is based on density rather than, 
as in our case, on the gravitational potential, as well as for 
haloes in the vicinity of clusters by JBasilakos ct al] (|2006l ) 
usi ng the moment of inertia e l lipsoi d of superclusters and 
by iRagone-Figueroa fc Plionij (120071 ) defining environment 
by the distance to the nearest cluster. It is clear from our 
present analysis that such alignments are in place at red- 
shifts of order one, and are maintained virtually unchanged 



performed with Ms/Mt = const., the filament direction is ob- 
tained on increasingly smaller comoving scales at higher redshifts. 
This partially erases the stronger correlation that is observed for 
the most massive haloes when the smoothing is kept at constant 
comoving scale for all redshifts. 



Figure 10. Median alignment angles between the halo major 
axis vectors and the eigenvectors pointing in the direction of fil- 
aments and normal to the sheets, respectively. Different redshifts 
are indicated with the three colours. Errorbars indicate the error 
in the median. The dotted line indicates the expectation value 
for a random signal. Data is shown for the ratio of the smoothing 
scale Ms/M, fixed. 

over the last eight or more billion years of evolution of struc- 
ture in the universe. 

3.5.2 Halo-Halo alignments 

We finally compute the spin-spin and spin-orbit correlation 
functions using the defini tions of IPorciani et all l|2002l ) and 
iBailin fc Steinmet j l|2005l ). While we show the results for the 
Ms/Mf = const, smoothing case, we stress that qualitatively 
the results remain unchanged when the constant smoothing 
is adopted. 

For the spin-spin correlation we have: 

aj(r) = {lJ(x).J(x + r)|), (12) 

where J is the intrinsic angular momentum of each halo, 
and the average is taken over all pairs of haloes which are 
separated by a distance r and reside in the same environment 
class. Similarly, the spin-orbit correlation is defined as: 

Cj.L(r) = {J(x).L(x + r)), (13) 

where L is the relative orbital angular momentum between 
two haloes separated by a distance r. 

Figure [11] shows the spin-spin alignment for haloes in 
clusters and filaments at the three redshifts of our study; up- 
per and lower panels show respectively the results for haloes 
with masses below and above M«. The shaded region shows 
the l(T-confidence area for the total sample, split by mass 
but not split by environment. The correlations within either 
of the environmental classes is never stronger than those 
for the total sample and all of them are consistent with no 
signal within 2a. Furthermore, we find no evidence for any 
significant redshift evolution of these correlations. 

The spin-orbit correlation function is shown in Figure 
1121 The strong correlation that we found at 2 = in Paper I, 
extending out to several Mpc, is present also out to redshift 
2 = 1 with no significant changes. 
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Figure 11. The mean alignment of intrinsic spin angular momentum between haloes in filaments and clusters at redshifts z = 1.05 
(left), 2 = 0.49 (centre) and 2 = (right). Data for filaments are shown in panels a) and c), clusters in panels b) and d). The upper 
panels a) and b) show the results for halo masses M < M*, and the lower panels c) and d) for halo masses M > M*. The dotted line 
indicates the expectation value for a random uncorrelated signal. The shaded region indicates the Icr confidence interval on the mean for 
the whole sample, split by mass but not split by environment. 
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Figure 12. The mean alignment of intrinsic spin and relative orbital angular momentum between haloes in filaments and clusters at 
redshifts z = 1.05 (left), z = 0.49 (centre) and 2 = (right). Data for filaments arc shown in panels a) and c), clusters in panels b) and 
d) . The upper panels a) and b) show the results for halo masses M < M« , and the lower panels c) and d) for halo masses M > M* . The 
dotted line indicates the expectation value for a random uncorrelated signal. The shaded region indicates the Icr confidence interval on 
the mean for the whole sample, split by mass but not split by environment. 



4 SUMMARY AND CONCLUSIONS 

We have used three N-body simulations, tailored to cover 
a range of almost five decades in mass with high resolu- 
tion haloes (> 300 particles), to investigate the dependence 
of halo shape, spin, formation redshift and alignment as a 
function of mass, environment and redshift. Using the tidal 
stability criterion of Paper 1 wo have classified haloes to re- 
side in four different cuvironments: clusters, filaments, sheets 
and voids. The attribution of haloes to these environments 
depends on one free parameter, Rg, the length scale used 
to smooth the underlying mass distribution. Relating this 
length scale to the mass contained in the Gaussian filter, 
Ms, in Paper 1 we optimised by visual inspection the red- 
shift zero value of = 10^=^ /r^M,, « 2M*(0), with M,(«) 
the typical mass scale collapsing gravitationally at redshift 
z. At the higher redshifts that wo study in this paper, wo dis- 
cuss two possible choices for the smoothing mass scale: i) a 
smoothing scale constant with redshift; and ii) a smoothing 
scale that varies such that Ms/M^ remains constant with 



redshift. The first approach leads to the median overden- 
sity in each environment increasing just as expected from 
non-lineax enhancement of density fluctuations; the second 
approach maintains the median density in each environment 
constant with redshift. 

In our analysis of the redshift evolution of the halo prop- 
erties we find that, when adopting a constant ratio M^/M* 
for the smoothing, the environmental influence is roughly 
invariant with redshift so that the mass scale at which the 
environmental influence sets in is roughly given by by the 
mass scale M,. Unveiling the importance of this physical 
mass scale in the onset of an environmental dependence of 
the halo properties is a first step towards understanding the 
origin of the environmental role in the evolution of dark 
matter haloes. 

Adopting the physically-motivated Ms/M^ = const, 
smoothing at all redshifts, we have investigated the depen- 
dence of the properties of isolated dark matter haloes of 
masses below and above the M, threshold on their environ- 
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ment. In general, we find that all halo properties show some 
dependence on environment for halo masses M < M, . The 
strength of the correlations, however, does not change much 
with redshift. There is virtually no redshift evolution of the 
halo properties when the correlations with halo mass are ex- 
pressed in terms of the normalised mass M/M* , indicating 
that the strongest evolution with redshift is related to the 
evolution of the mass scale for collapse, M, . In detail, our 
main results are summarised as follows: 

• There is a strong environmental dependence of halo 
formation times with environment for haloes with masses 
M < M«. At any given mass in this mass regime, haloes in 
clusters tend to be older than haloes in the other environ- 
ments, and haloes in voids form much more recently than in 
any other environment. The global halo population (with no 
splitting for environment) is well described by a logarithmic 
relation between mass and formation redshift with a roughly 
constant slope with redshift. 

• The median spin parameter X[^c^ of the total halo sam- 
ple, not split by environment, increases weakly with mass up 
to around lOA/* at all redshifts. There is no significant resid- 
ual dependence on redshift besides the mass rescaling with 
M«. There is an additional tendency for high mass haloes 
in filaments to spin more rapidly than haloes of the same 
mass in clusters. Haloes in voids have the lowest median 
spin parameters. 

• Haloes of a given ratio A//A/* have very similar median 
shape parameters independent of redshift in the mass and 
redshift ranges that we have investigated. Independent of 
environment, haloes are increasingly more spherical and less 
triaxial the lower their mass. Haloes with masses M < Af* 
are slightly more spherical and more oblate in clusters than 
in filaments, and there is a hint that the situation reverses 
for M > Mt, i.e., that haloes in filaments are more oblate 
than cluster halos at high masses. Low-mass haloes in voids 
have systematically the lowest median sphericity of similar 
mass halos in denser environments. 

• In the M < Af* mass regime, haloes in sheets tend to 
have spin vectors in the plane of the sheets, and haloes in 
filaments tend to have spin vectors pointing along the fil- 
aments; above the Af, mass scale, there is evidence that 
haloes in filaments reverse the previous trend and tend to 
have spins orthogonal to the filaments. Furthermore, haloes 
with masses M > O.IA^* tend to have their major axis paral- 
lel to their host sheets or filaments, with the strength of the 
alignment increasing with increasing mass. This may reflect 
the fact that, for massive haloes, the gravitational poten- 
tial field is substantially infiuenced by their shape and thus 
leads to an aligned tidal field. The alignment of halo spins 
and major axes with the large-scale structure that we have 
unveiled up to redshifts of order z = 1 should be taken into 
account in studies o f weak lensing maps of cosmic shear (eg. 
ICatelan et aLlfioOll ). especially in sheets and thus in regions 
surrounding voids. 

• There is no evidence for a significant spin-spin corre- 
lation between neighbouring haloes. There is in contrast a 
substantial halo spin-orbit alignment, whose strength ap- 
pears to be independent of mass, environment and redshift 
up to z ~ 1: haloes in close pairs tend to spin preferentially 
parallel to the orbital angular momentum of the pair. 

An important conclusion that we draw from our study is 



that the environmental infiuence on halo properties shows 
an intriguing dependence on the halo mass, and appears to 
be essentially modulated by the typical mass scale of gravi- 
tational collapse Af* at each redshift. Our data suggests that 
the Af = Af, (z) mass scale might indeed play the role of a 
bifurcation point below which many of the median proper- 
ties of dark matter haloes either begin to feel the influence of 
their large-scale environment, or show an opposite response 
to their large-scale environment relative to the more massive 
haloes. The existence of such a thresholding mass scale in 
the environment-halo relationship is yet to be understood. 
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